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Abstract 

We develop five new statistical measures of the quality of fits, which we 
combine with the usual confidence level to determine the models which fit 
best all available data for total cross sections and for the real part of the 
forward hadronic amplitude. 



"'^presented by J.R. Cudell at the 6*'* workshop on non-pcrturbative QCD, American 
University of Paris, 5-9 June 2001. 
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Phenomenological studies compare models with data, in order to deter- 
mine whether these agree together. The typical criterion is that there is good 
agreement if the per degree of freedom {x^/dof) is of the order of 1, or if 
the confidence level is bigger than a typical value. To take an explicit exam- 
ple, one can imagine fitting all available elastic hadron-hadron amplitudes 
(their imaginary part being provided by the total cross sections, and their 
real part by the p parameters) to analytic parametrisations|]l[. Although we 
shall only treat this case explicitly, the remarks and tools given here do apply 
to most situations. 

There are several problems with this conventional approach: first of all, 
models usually do not apply everywhere, hence the comparison should hold 
only for part of the data. In the case of hadronic amplitudes, the models 
are smooth analytic functions, which are expected to work at high enough 
energy. Some of these models are also expected to have large (unitarising) 
corrections at large energies. Hence in this simple case, there is a range of 
energy over which the comparison should be meaningful. However, as is often 
the case, the exact range over which the models hold is not predicted by the 
theory, and should emerge from the fits themselves. 

Secondly, once one fits many points, a large discrepancy between the 
theory and a few of the data points can be overshadowed by a good overall 
agreement. This has the drawback that it is precisely these points which may 
point out to new physics, but at the same time this may reveal problems with 
the data. Hence some uniformity in the description of the data is needed, 
and studies such as the present one, applied to soft hadronic amplitudes, lead 
to a reassessment of the data used in the fit. 

Thirdly, the fits to models lead to some values for the physical parameters 
of the model. Here, one must take into account the interplay between these 
values and the data sub-sample for which the model applies. This means 
that if one changes that sample slightly, the parameters extracted for the 
models should be stable. A typical counter-example is found in the fits to 
a simple-pole pomeron, which give (wrongly) an increasing pomeron intercept 
once the minimum energy of the data is below 9 GeV. 

Because we want to be able to consider a large data sample and many 
possible models, we are aiming at the development of an automatic decision- 
making procedure, and hence we want to measure the above criteria. Al- 
though we are not entirely finished with this program, we can present what 
seems a reasonable set of measures which reflect the above aspects of the 
fits. All these measures, or indicators, are constructed so that the higher 
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their value the better is the quahty of the data description. 

The first indicator concerns the sample of data that can be fitted. In the 
case of hadronic amplitudes, we shall consider the range of energies where 
a given model applies, defined as the region in energy where the fit has a 
confidence level (CL) bigger than 50%. Its size will be one of the measures 
of the quality of the fit: we define the applicability A of model M as: 



nM,high\ J2 ' A^ 

^M,low - ^^^^ 



= -.log|^i^h ^"-^ (1) 

where E^''^^^^ ( resp. E^'^°^ ) are respectively the highest and lowest vahics 
of the energy in the area of applicability of model M in the data subset j 
and Wj is the weight determined from the best fit in the same interval. After 
we have defined where the model may work, we can check how well it fits, 
although by definition all models will provide a satisfactory fit. We may 
consider the usual confidence level, Cf^ = CL{%), where the CL refers to 
the whole area of applicability of the model M, or a reduced one limited 
to the intersection of the areas of applicability of all models qualified for the 
comparison. 

The next measure of quality has to do with the number of parameters 
of the model, given the number of data points in the range of applicability. 
Hence we define the rigidity Ri as: 

' par 

All the information on a given fit is contained in the error matrix, and we 
shall use it to define the new measures. Hence the first condition to check 
is whether the error matrix itself is reliable, i.e. whether the correlations 
between parameters are minimal. Hence we define the reliability R2 as: 

R¥ - ^ ' • E 0(90-0 - c^) (3) 

where C/lf is the correlation matrix element in % calculated in the fit at the 
low edge of the applicability area. 

We are now in a position to define the stability of the model with respect 
to the data range considered. In the case of hadronic amplitudes, three pos- 
sible changes can be considered: we can consider the variation that comes 
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from modifying the energy threshold of the fit (energy stability Si), or 
the fluctuation of the from bin to bin for some data binning motivated 
by physics (in the case of hadronic amplitudes, we bin according to the 
process/observable) (uniformity U), or the reproducibility of the param- 
eters values when fitting, with the same number of adjustable parameters, 
a reduced data sample and a reduced number of observables, in the case 
of hadronic amplitudes when excluding the real part (r-stability S2)- The 
latter is introduced in this case because the data for p parameter data may 
be less reliable than those for the cross section. Hence we obtain the three 
measures: 
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where: P* is the vector of parameters values obtained from the model fit 
to the whole area of applicability; P**'^?' is the vector of parameters values 
obtained from the model fit to the reduced data set on the step, in our case 
step means a shift in the low edge of the fit interval to the right by 1 GeV; 

and are the error matrix estimates obtained from the fits to the 

total and to the reduced data samples from the domain of applicability, t 
denotes the total area of applicability, and t{no p) the data sample with p 
data excluded. 

Having defined these measures of the quality of fits, we want to use them 
to see whether we can decide which is the safest model to use to reproduce 
a given set of data. As already emphasized in the case of hadronic ampli- 
tudes, standard methods do not allow one to decide which models are to be 
preferred, and several classes of parametrisations are possible. Using these 
new measures, we can try to decide which models are best. All models con- 
sidered are the sum of several terms: the low-energy sector is described by 
an amplitude with the following imaginary part, with Si = 1 GeV^: 

/m(A"'') = Y^'' {s/siY' T (s/si)"' (5) 

The first term has charge-conjugation C = +1 whereas the second has 
C = —1 (with the — sign for a positively charged beam). These two terms are 
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symbolized by the notation RR in the following. The high-energy behaviour 
is dominated by a pomeron term, for which we consider the following terms, 
or their combination, in the imaginary part of the amplitude: 



ImiA"'') 
Im{A''^) 



(6) 
(7) 
(8) 
(9) 



= S»^sln(s/si) 



which we denote respectively hy E, P, L and L2. If we take for the pomeron 
a simple pole model we obtain in this notation RRE, whereas a double pole 
gives RRPL and a triple pole RRPLL2, which we can write as RRPL2 with 



Furthermore, we have considered several possibilities to constrain the 
parameters. The following notations are attached as either superscript or 
subscript to the model variants in each case: 
d means degenerate leading reggeon trajectories ai — 
u means universal (independent of projectile hadron); 

nf means that we have not imposed factorization for the residues of the 
pomeron term(s) in the case of the 77 and jp cross sections; 
qc means that a quark counting rule is imposed on the residues of the am- 
plitude for Ep scattering, constrained by the residues in pp and Kp; 
c implies the use of the Johnson- Treiman-Freund relation for the cross sec- 
tion differences: Ao-(A^) = 5A(T{n), Acr{K) = 2Ao-(7r). 

Finally, the real parts of the amplitudes can be obtained through s —>■ u 
crossing. 

The first results are quite generic and based on a study of the alone: 

(1) All analytic descriptions of the data based on the above terms break 
down at i/s < 4 GeV; 

(2) Most models require a non-degeneracy of lower trajectories. Degeneracy 

can be accommodated only by RRPL2; 

(3) Simple pole pomerons fail to reproduce the real part of the cross section, 

and all models have problems in reproducing some of the p data; 

(4) Cosmic ray data are well reproduced by the best parametrisation, with 

no need of re-analysis of the published data; 



a scale Sq instead of Si in the log^. 
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Table 1: Best models for total cross sections 





ACCURRS 


ACCURRS20 


AURS20 


1 




RqcRcL'^qc 


Rqc^c^^qc 


2 


{RR)dPL2u 


{RRc)dPL2u 


{RRc)dPqcL2u 


3 


{RRc)dPL2u 


{RR)dPL2^ 


RqcRcL qc 


4 


Rqc-^c-^'^qc 


RRL2qc 


{RRc)dPL2u 



(5) Although quark counting rules can be approximately implemented, it is 
also possible to have a universal rising term for the pomeronP|. 

Furthermore, we can use all the information contained in our indicators to 
define the best models. Several schemes are possible: 

i) the ACCURRSS scheme: we take all indicators, including the Ci and C2, 
and for each indicator we order the models considered from rank 1 to A^ 
according to the value of the indicator. We then sum the 8 numbers obtained, 
and the best model is the one with highest rank overall; 

ii) As the indicators are statistical measures, we can do the same as above 
but consider that a model is better than another (and give it one point) only 
if its indicator is bigger than that of the other model by e.g. 20%. This leads 
to the ACCURRSS20 scheme; 

iii) Finally, one may argue that all the CL are acceptable, and that the 
number of parameters is not relevant as we have chosen functional forms, 
and hence in principle an infinite number of parameters. Using a statistical 
ranking similar to ii) leads then to the AURSS20 scheme. 

Using these ranking schemes, we obtain the best models (out of the order 
of 30 variations on the terms used in the models) given in Table 1 for fits 
to total cross sections only, and in Table 2 for fits to all data for hadronic 
amplitudes. As can be seen from these table, simple-pole pomerons are never 
preferred, and models containing a log^ s rise in the cross section always pro- 
vide the best fits to the data. To reach a more restrictive conclusion, one 
needs to use the p data, in which case the preferred model is consistently 
RRPL2u. The problem however is that the p data are poorly reproduced 
by all models considered in this study, hence one cannot be sure that this 
preference will survive future iterations of the cross assessments with new 
models and new data added, such as off-forward cross sections, or DIS struc- 
ture functions, or data from future experiments. As is always the case in 
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Table 2: Best models for hadronic amplitudes 





ACCURRSS 


ACCURRSS20 


AURSS20 


1 


RRPL2u 


RRPL2u 


RRPL2u 


2 


RRLnf 


RRcL2qc = RRL2qc 


RRL2 


3 


{RRc)dPL2u 




4 


{RR)dPL2u 


{RR)dPL2u = RqcRcL2qc = RRL2 


{RR)dPL2u 



these studies, it would be of utmost interest to have higher-energy data for 
other beams than protons and antiprotons. This would of course enable one 
to determine directly whether the pomeron counts quarks, or has a universal 
component. 
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